Persistent Current in the Ferromagnetic Kondo Lattice Model 
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In this paper, we study the zero temperature persistent current in a ferromagnetic Kondo lattice 
model in the strong coupling limit. In this model, there are spontaneous spin textures at some values 
of the external magnetic flux. These spin textures contribute a geometric flux, which can induce 
an additional spontaneous persistent current. Since this spin texture changes with the external 
magnetic flux, we find that there is an anomalous persistent current in some region of magnetic flux: 
near <E > /<I > o = for an even number of electrons and $/$() = 1/2 for an odd number of electrons. 



I. INTRODUCTION AND MODEL 

The Berry phase! plays important roles in transport 
in mesoscopic systems. The most well known effect is 
the persistent current in quasi-one-dimensional nor 
metal mesoscopic rings threaded with magnetic flu: 
This purely quantum effect manifests itself through the 
acquisition of a phase factor proportional to the magnetic 
flux which can change the boundary condition of the or- 
bital wavefunctionax Similarly, when a quantum spin 
adiabatically follows a magnetic field varying in space or 
time, the phase of its wavefunctions acquires an addi- 
tional contribution - a geometric or Berry phase. It was 
shown by Loss et al. that this geometric-phase can also 
induce persistent spin and mass currents! 

In this paper, we study the persistent current in the 
ID ferromagnetic Kondo lattice model (FKL) in the 
strong coupling limit. The FKL model in the strong co 
pling limit is related to the double exchange mechanisn 
which has been studied recentlyuii3 due toits relevance 
to the colossal magnetoresistive Mn-oxidesLj. From now 
on, we will call the FKL model in the strong coupling 
limit the double exchange model. In this double exchange 
model, the ground state has spin textures at some val- 
ues of magnetic fluxtrta Since the electronic spins are 
strongly ferromagnetically coupled to the local spins, the 
spin texture of the local spin can induce a Berry phase 
in the electronic wavefunctions. The geometric flux pro- 
duced by the local spin texture then induces spontaneous 
persistent currents. In previousljj-,studied persistent cur- 
rents induced by geometric fluxO, the geometric Berry 
phase is due to adiabatically varying external magnetic 
field textures. In the double exchange model, the geomet- 
ric Berry phase is due to the local spin textures, which 
depend self-consistently on the dynamics of the electrons. 
With the change of external magnetic flux, the local spin 
texture varies. This in turn will change the geometric 
phase-induced persistent current. 

Here we study the zero temperature persistent current 
of an ideal one-channel double exchange ring threaded 
with a magnetic flux <F The effects of finite temperature 
and disorder (elastic scattering) for double exchange sys- 



tems are similar to those of the conventional persistent 
current in normal metal ringstJ, and will not be consid- 
ered here. The FKL model for a A-site ID ring with 
magnetic flux $ can be written as 



h.c. 
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where <I>o = hc/e, a is the Pauli matrix, and Si is the local 
spin. The operators c ia (cj CT ) annihilate (create) a mobile 
electron with spin a. For the double exchange model, 
Ju/t ^rjl- The zero temperature persistent current is 
equal toQ I — —cdE/d<p. It is easy to show that the 
free energy and persistent current are periodic in <I> / <!>□ . 
In Hamiltonian ([!]), we have neglected electron-electron 
interactions. In the double exchange model, the on-site 
electron-electron interaction is unimportant; it mainly 
renormalizes Jh- 

In Sec. |fi], we will study the properties of spin tex- 
tures and persistent current in the classical spin limit. In 
this limit, many results can be obtained analytically. In 
Sec. [II, we study the effects of quantum fluctuations of 
order 1/S using a set of trial wavefunctions. In Sec. IV, 
we study the spiral ordering and persistent current for 
5=1/2 rings using exact diagonalization. We summa- 
rize our results in Sec. 
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FIG. f. Local spin configuration with spiral ordering 
nj = {OjAj) = (8, 2 <!*>)■ 
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FIG. 2. Two equivalent local spin configurations. 



II. CLASSICAL SPIN LIMIT S -> oo 

The zero temperature persistent current is most eas- 
ily calculable in the classical local spin limit. In this 
limit, the electrons are moving in a frozen local spin 
background. Hamiltonian ([!]) is then transformed to a 
one-body Hamiltonian. The ground state of Hamiltonian 
(|l|) is the ground state with optimal static local spin con- 
figuration. Since the original Hamiltonian is translation 
invariant and there is an effective short range ferromag- 
netic coupling between local spins due to the double ex- 
change mechanismfl, we conjecture that the ground state 
has spiral ordering 



hj = (9j,<l>j) = (9,j<t>o), 



(2) 



where <fio = 2nnr/N is an integral multiple of 2ir/N [see 
Fig. [TJ . Note that the maximally polarized ferromagnetic 
ordering is an extreme case of this "spiral ordering" with 
6 = 0. Here 9 is a variational parameter. We need to 
minimize the total energy E to t{9) to find the optimal 9 
for the local spin configurations. 

Using the local spin configurations of Eq.(||), the 
Hamiltonian in fc-space is 
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where Vs = A k sin6>, A k = tcos(k — O — w§-)~ tcos(k- 

T x = e° - A fc cos0 and T 2 = e° 
the kinetic energies with e k = t cos(fc 



z k t A k cos 9 are 
M;) + tcos(k~ 

4>o — From Eq.@, we can show that the electronic 

energies with the two spin configurations 9 <-> n — 9 1 
cj)Q <-* — 0o are degenerate [see Fig. 



The lower eigenenergy of the single electron states are 



SJh 



y 2 + A 2 k -A k SJ H cos9. 



(4) 



For SJ H > |A 2 |, 
SJh 



Ck 



e° k + A k cos 9 ~ A 2 k sin 2 9/ (SJ H )- (5) 



Note that the eigenenergy is periodic in $/$()■ 

In the limit S = oo, the many-body wavefunction is 
just the Slater determinant of single particle wavefunc- 
tions. The total energy is 



E, 



(6) 
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At zero flux <E>/<I>o = 0, the ground state fat an even num- 
ber of electrons has spiral spin orderingO 9 = w/2, i.e. 
the spins lie on the equator with zero total magnetization. 
For densities not close to half filling n = N (n/N < 0.8 
for SJh ~ 124), the lowest energy spiral state has a sin- 
gle twist: 0o = 2tt/N. Actually, with zero flux, all the 
spiral states with finite 9 and a single twist have lower 
energy than the lowest ferromagnetic state (9 = 0). The 
energy difference is: 



E to t(9) - E tot (0) ~ - sin 2 9 /(SJh) 



\k\<k f 
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for SJh ^S> |A||. The results for an odd number of elec- 
trons is similar but shifted by $/$o — * $/$o + l/2. Here 
the Berry phase plays an important role to lower the en- 
ergy of spiral states: with the spiral spin texture, the 
magnitude of the effective hopping is reduced, but the 
effective hopping matrix element has an additional com- 
plex phase, which can lower the total energy of spiral 
states for even (odd) number of electrons at <&/<!> o = 

(*/#o = 1/2). 

Let us consider what happens if there is an infinitcsi- 
For SJ H > 



mal external flux f - < 1 



IA 2 fe | 
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E t ot(9) - E tot (0) ~ 2(1 - a»e)twak f sin( — — ) 
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Thus for small external flux ^- < 1, the spiral angle 9 
decreases as 



cos l 



oc sin( 
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and we can expect that as the magnetic flux &/$o 
changes, the pitch of the spiral ordering remains constant 
(i.e. 0o = 27r/iV), and the angle 9 of the spiral ordering 
decreases continuously. At some critical value $ = $ c , 
9 = 0, i.e. the spins become ferromagnetically ordered. 
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In Fig. |2|, we show the change of spiral ordering angle 9 
with magnetic flux. 
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FIG. 5. Density dependence of the critical magnetic flux 
<J> C in the classical spin limit. SJh = 12£, N = 300. 



FIG. 3. Change of the spiral ordering angle 6 with magnetic 
flux $/$ at SJh = 12t and 5 = oo for a N = 450 ring with 
n = 300. 

At fixed density, the critical value of magnetic flux $ c 
depends only weakly on the size of the ring N, as shown 
in Fig. [|. From Fig. [|, we can see that the critical value 
increases slightly with increasing system size N and satu- 
rates at N < 10 3 . As shown in Fig. |[ $ c depends nearly 
quadratically on the electron density. From this density 
dependence, we can also see that there is no particle-hole 
symmetry in this model. Since at half filling n — N the 
ground state is antiferromagnetically ordered, our calcu- 
lations are not valid close to half-filling. We also show the 
SJh dependence of <& c in Fig. ^, from which we can con- 
clude that it is very close to a power law for large SJh- 
$ c ~ C/(SJh)- This 1/(5 J/f) dependence is expected 
from Eq. (||). 
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FIG. 6. Hund's rule coupling dependence of the critical 
magnetic flux $ c in the classical spin limit, n = 200, N = 300. 
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FIG. 4. System size dependence of the critical magnetic 
flux $ c in the classical spin limit. SJh = 12t, 24t, density 
v = 2/3. 



When the local spins are not ferromagnetically or- 
dered, their spin texture will induce a Berry phase in 
the electron wavef unctions. Thus finite 9 will induce a 
geometric flux, which in turn induces a spontaneous per- 
sistent current in the mesoscopic double exchange ring. 
This spontaneous persistent current will cancel part of 
the persistent current due to finite external magnetic 
flux. In Fig. (j^) , we show structures of the persistent 
current for even and odd numbers of electrons. In con- 
ventional mesoscopic systems, the zero temperature per- 
sistent current is linear in <I> for fc$o < $ < (fc + l)$o 
for an even number of electrons, and (k — l/2)$o < $ < 
(k + l/2)$o for an odd number of electrons. In a meso- 
scopic double exchange ring, due to the additional contri- 
bution to the persistent current from spin textures, the 
persistent currents are no longer monotonic in the regions 
mentioned above. Upon ensemble averaging of rings with 
different electron numbers, the period of the persistent 
current is changed from $ = $o to $ = $o/2. For dou- 
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ble exchange rings, the persistent current has anomalous 
regions near both $ = $o and $ = <&o/2 after ensem- 
ble averaging. In the remainder of this paper, we will 
continue to discuss even and odd numbers of electrons 
separately. 
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FIG. 7. Persistent current in a double exchange ring in 
the classical spin limit: N = 450, n — 300 (solid line), and 
n = 301 (dashed line), SJh = 16. The dotted-dashed lines are 
the persistent currents in a conventional mesoscopic metallic 
ring. 



III. SEMICLASSICAL SPIN LIMIT S > 1 

If the local spin S is large but finite, the local spins 
have quantum fluctuations of order 1/S 1 , which can be 
taken advantage of by electrons to lower the electronic 
energy. This is the spin polaropic effect, which is similar 
to that in the Hubbard modelia. Here we still consider 
the case where the ground state is either spirally ordered 
or ferromagnetically ordered, defined by hj = (9j ,<j>j) — 
(9,j(j)o), as in the classical limit. 

First we consider the eigenstates of the single itinerant 
electron case. If one electron is at site j, the complete 
spin basis for this site can be chosen as 
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with a quantization axis of hj. In the semiclassical limit, 
these total spin states \S ± ^, S ± ^ — k) with the local 
quantization axis of hj are an appropriate basis to de- 
scribe the spin states at the site with the itinerant elec- 
tron present. In the classical limit as studied in the previ- 
ous section, only \S± \, S± ^) need be used to construct 
the exact wavef unctions. In the semiclassical limit, the 
quantum fluctuation corrections will come from the con- 
tributions of \S ± |, S ± | - k) (k > 0). If we include 



just the first order quantum fluctuation effects, we need 
a four trial wavefunction basis: 



h#>> 
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where | • • •) - is the spin state | • • •) with the local quan- 

n ¥ i . . . 
tization axis hj , which is a generalized coherent stateO in 

the global quantization axis. Although we cannot show 
that the trial wavefunction basis in Eq. ( jll| ) rigorously de- 
scribes the 1/S corrections to the results in the classical 
spin limit, we believe this will at least describe the qual- 
itative effects of 1/S corrections. Note also that one can 
include a variational component in the empty sites wave- 
functions such as { | S, S) 



riift 



c\S,S-l} ~ }, where c 



is a variational parameter. But one can show that these 
corrections to the eigenenergy will be higher order in 1/S* 
and thus negligible. 

Using the variational basis (|ll|), the Hamiltonian in k- 
space can be readily written down. For simplicity, we will 
only discuss the large Hund's rule coupling limit SJh = 
oo. Then we only need the basis states \ipj ) and IV^)- 
The Hamiltonian in this variational 2x2 basis is 



H = 




kiVs 



2S+ 



SJh 
2 



(12) 



where T 3 = ^i+l( e k + A fc cos(9). 

From the above trial wavefunctions we can see that 
quantum fluctuations make the spins near the electron 
more polarized than the average spiral ordering direc- 
tion. Thus we can think of the quantum fluctuation ef- 
fect as a spin polaron effect. In the limit Jh — > oo, the 
eigenenergies of the single electron states are 



SJ 



H 



V 



A fc cos 9 + 5) 



~yj (e° A k cos 6 6)* + ^^^t sin 2 9, (13) 

where 5 = 1 ^ TI {e° k +A k cos9). When |e£-A fe cos 9-6\ > 
\2A k sin0/|V2S + 1, 

SJh o a a A 2 k sm 2 9 
ek = -—- ek + AkCOSd - 2S\e°-A k cos9-6\ - 

(14) 

In the semiclassical limit, we can construct approx- 
imate many-body wavefunctions using Slater determi- 
nants of single particle wavefunctions. However, this 
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construction is only approximate because we neglect the 
interaction between different spin polarons. The correc- 
tions due to the spin polaron interaction is 0(n/S 2 ). 
Thus this approximation is only valid for the dilute elec- 
tron limit and large S. Similar to the finite Jh effects, 
the quantum fluctuations induce a spiral ordering in- 
stability at $/$o = for an even number of electrons 
and $/$o = 1/2 for an odd number of electrons. From 
Eq. (pTh , the lowering of the GS energy with spiral order- 
ing due to quantum fluctuations is E sp — Efm oc — \t\/S. 
Figure (g) show the persistent current for semiclassical 
spins. For low electron densities, the 1/ S quantum fluc- 
tuation effects are similar to those of finite SJh in the 
classical spin limit. 
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FIG. 8. Persistent current in a double exchange ring in the 
semiclassical spin limit: N = 300, n = 120 (solid line), and 
n = 121 (dashed line), S = 5, Jh = oo. The grey long-dashed 
lines are the persistent currents in a conventional mesoscopic 
metallic ring. 



IV. PERSISTENT CURRENT IN 5= 1/2 
SYSTEMS -EXACT DIAGONALIZATION STUDY 

From the results in the classical and semiclassical lim- 
its, we expect that the spiral ordering instability will be 
enhanced in the quantum spin limit. This instability for 
general spin S has been demonstrated previouslytij us- 
ing a many-body trial wavefunction. To confirm that the 
ground state spin texture is indeed spiral for a quantum 
spin, we calculate the exact ground state for S = 1/2 
rings using exact diagonalizationlij. For an even num- 
ber of electrons, we find that the ground state has total 
spin Stotai = for both large SJh and SJh = oo at 
$/$o = 0- With an increase of external flux, the to- 
tal spin Stotai of the ground state increases, to Stotai = 
gmax = ( NS + n / 2 ) a t $/$ = 1/2. For an odd 
number of electrons, the results are shifted by $/$o — ¥ 
$/$o + 1/2. We speculate that the non-ferromagnetic 



state with S to tai < S max is a sum over global rotations 
of spiral states with polarization angle 9 (0 < 9 < tt/2). 
To check this statement numerically, we calculate the cor- 
relation function ((Si x S2) • (Si x Si+i)) for even numbers 
of electrons at $/$o = 0. Because this ground state with 
Stotai = can be thought of as a sum over all global ro- 
tations, we need to look at scalar correlation functions. 
The most relevant correlation function for the spiral or- 
dering is ((Si x 5 2 ) • (Si x Si+i)). In Fig.(9R, we show 
the spiral correlation function ((Si x S2) ■ (Si x Sj+i)) 
for different size rings of spin S = 1/2. The corre- 
lation is only weakly dependent on i and has magni- 
tude ~ (0.25 sin(27r/iV)) 2 . This correlation function for 
S = 1/2 is closer to (0.25 sin(27r/7V)) 2 as the system 
size N and/or electron number n increases. This clearly 
demonstrates that the spin ordering is spiral, similar to 
what we find in the classical and semiclassical limits. 

The geometrical flux due to the local spin textures also 
contributes to the persistent currents in quantum spin 
systems. In Fig. (|l(i|), we show the exact diagonalization 
results for the persistent current in S = 1/2 double ex- 
change rings. We can see from Fig. ( |Io|) that the width 
of the anomalous region is increased compared to that 
in the classical [see Fig. (0)] and semiclassical limit [see 
Fig. (||)]. This demonstrates that the possibility of spi- 
ral spin ordering is enhanced in the quantum spin limit, 
as expected from the 1/5 results. Although we do not 
have enough numerical results for $ c to study the depen- 
dences on system size N, electron density n, and SJh, 
we believe that these dependences will be similar to those 
for classical local spin double exchange rings, as shown 
in Figs.(@)-(|). 
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FIG. 9. Correlation ((Si x S 2 ) ■ (Si x 5<+i)) for S = 1/2 
rings normalized by (0.25 sin(27r/iV)) 2 . SJh = 40t. 
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FIG. 10. Persistent currents in the 5=1/2 double ex- 
change system. N=6,8 and n=2,3. SJh = 40t. 



V. CONCLUSION 

In conclusion, we have studied the novel zero temper- 
ature persistent current in the double exchange system. 
In this system, there are spontaneous spin textures at 
some values of the external magnetic flux. These spin 
textures contribute a geometric flux, which can induce 
additional spontaneous persistent current. Since the spin 
textures vary with the change of external magnetic flux, 
there are anomalous persistent currents in the region 
near <&/<!> o — for an even number of electrons and 
$/3>o — 1/2 for an odd number of electrons. After en- 
semble averaging of double exchange rings with different 
electron numbers, the persistent current has a period of 
$ = ^o/2, and the anomaly occurs at 3>/<I>o = and 1/2. 
The effect of disorder (elastic scattering) and finite tem- 
perature for the double exchange systems are similar to 
those in conventional normal metal ringsEj. With the ad- 
vance of fabrication technology, this fascinating anoma- 
lous persistent current and the magnetic flux dependent 
spin texture may become observable. 
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